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a comparison to external galaxies 
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Abstract The Milky Way bulge offers a unique opportunity to investigate in detail 
the role that different processes such as dynamical instabilities, hierarchical merg¬ 
ing, and dissipational collapse may have played in the history of the Galaxy forma¬ 
tion and evolution based on its resolved stellar population properties. Large obser¬ 
vation programmes and surveys of the bulge are providing for the first time a look 
into the global view of the Milky Way bulge that can be compared with the bulges of 
other galaxies, and be used as a template for detailed comparison with models. The 
Milky Way has been shown to have a box/peanut (B/P) bulge and recent evidence 
seems to suggest the presence of an additional spheroidal component. In this review 
we summarise the global chemical abundances, kinematics and structural properties 
that allow us to disentangle these multiple components and provide constraints to 
understand their origin. The investigation of both detailed and global properties of 
the bulge now provide us with the opportunity to characterise the bulge as observed 
in models, and to place the mixed component bulge scenario in the general con¬ 
text of external galaxies. When writing this review, we considered the perspectives 
of researchers working with the Milky Way and researchers working with external 
galaxies. It is an attempt to approach both communities for a fruitful exchange of 
ideas. 
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1 Introduction 


What is the origin of the bulge of the Milky Way? The answer to this question is 
a crucial step towards identifying the history of events that took place during the 
formation and evolution of the Galaxy as a whole. As a matter of fact, the answer 
to this question has changed dramatically, since the early times of Galactic bulge 
archaeology until today, between the two main ideas behind bulge formation: the 
merger-driven bulge scenario, where a bulge is formed violently and quickly during 
the early stages of the Galaxy dominated by the gravitational collapse or hierar¬ 
chical merging of sub-clumps of dark matter carrying baryons and gas (e.g. |Abadi| 
et al 2003 1 Elmegreeri| |1999| >, and the secular evolution scenario where the bulge 
structure is naturally born from the dynamical evolution of the stellar Galactic disc 


( Combes and Sanders| |198l] |Raha et al[|199rt|Norman et al| |T996[ |Athanassoula| 

2005] >. Until this day, we cannot say with absolute certainty which one of these sce¬ 
narios, or ifperhaps both, played a major role during the formation of the bulge of 
our Galaxjrl 

The reason for this long term debate might lie, ironically, in our greatest strength. 
The Galactic bulge allows us to investigate its properties by taking advantage of the 
fully resolved stellar populations - a unique strength that can be understood as the 
ability to see its properties in a unique level of detail with respect to what we can 
learn from the observation of external bulges. However, such an advantage also 
means that nearly 500 sq. deg. of sky must be homogeneously covered in order to 
obtain the most global picture of the Bulge. To obtain a general characterisation of 
these properties is a crucial step in order to answer the question of the origin of our 
Galactic bulge. 

As a consequence, observational efforts during the last decade have been focused 
on solving this limitation. As a result, our knowledge regarding the global properties 
of the Galactic bulge has increased considerably thanks to the advent of dedicated 
spectroscopic and photometric surveys. We are currently witnessing a revolution in 
the field of Milky Way bulge research that will also find its place within the bulges 
of other disc galaxies. It is thus a moment in which the communities of Galactic 
and Extragalactic research are approaching each other. From this, one cannot expect 
anything but a fruitful exchange of ideas that will certainly push both fields forward. 
However, it is not straightforward for members of each community to study the other 
held, for at least two reasons. The first obstacle is the vast amount of work, rich in 


1 To be complete, we must mention the pioneering work of Eggen et al 1962), who suggested that 
the first stage in the formation of the Milky Way was a fast (< 10 s yr) monolithic collapse of a 
single massive gas cloud, which could have formed both the Galaxy stellar halo and the bulge. This 
scenario was later replaced with the merger-driven scenario mostly due to the widely varying ages 
of different components of the Galaxy and the hierarchical nature of LCDM theory. We refer the 
reader to Brooks & Christensen (this volume) for a review on merger-driven bulge formation. In the 
last decade or two, the secular evolution scenario has slowly, but unequivocally, gained terrain over 
the merger-driven scenario. Another bulge building scenario we do not discuss here has recently 
been put forward by e.g. |Elmegreen etal|j2008| - see review by Gadotti 2012). In this scenario, 
bulges form by the coalescence of giant clumps in primordial discs. This scenario can explain the 
formation of spheroids but does not account for box/peanuts. 
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details, that one has to become familiarised with. The second obstacle is the jargon 
employed independently by each group which hampers understanding. This review 
is a modest first attempt to overcome these obstacles. 


2 The structure of the Milky Way bulge 


The most basic definition of a galactic bulge is that of an over-density that swells up 
from the plane of the disc. The idea of this natural conception originated from the 
observations of other disc galaxies, in particular the bulges of edge-on spirals, which 
allow us to compare them more easily with our own in a purely morphological way. 

Within spiral galaxies, we could in principle distinguish between classical bulges 
and those bulges which are formed via secular evolution, based almost solely on 
their morphological signatures. Clearly, a proper characterisation of the structural 
properties of the Milky Way bulge would provide us with a valuable set of con¬ 
straints needed to find its place in the more general scheme of external bulges and, 
as extensively shown by galaxy formation models, to connect these constraints to 
the different mechanisms of origin. 

After it was first postulated by de Vaucouleurs (1964j>, followed by Sinha (1979|> 
and Liszt and Burton (1980) among others, Blitz and Spergel (1991) predicted the 
presence of the bar-like structure for the inner regions of the Galaxy based on in¬ 
frared observations. Later on, continuing with the exploitation of infrared imag¬ 
ing in order to overcome the strong dust obscuration towards the inner galaxy, the 
COBE/Diffuse Infrared Background Experiment ( [Smith et~al 2004| > data was used 
by Weiland et al (1994 1 to unambiguously establish the presence of the bar. Soon 
after, the COBE data further revealed the global B/P morphology of the Milky Way 
bulge (Dwek et al, |1995) >. 

An important number of Bulge structural studies have been based on the stel¬ 
lar counts of red-clump stars, which are the metal-rich counterpart of the well 
known globular cluster horizontal-branch stars. The absolute magnitudes of red- 
clump stars, are found to have little dependence on age and metallicity, making them 
one of the most powerful tools for deriving distances towards the bulge and there¬ 
fore tracing its global morphology. This method is based on the construction of the 
luminosity function of the Bulge towards a given line of sight where the red-clump 
feature can be easily identified and fitted with a Gaussian distribution to obtain the 
mean red-clump magnitude (Stanek et al |1994| l. Zoccali ( 201 0[ > and |Mc W illiam 
|et al| ( |2010| l presented the discovery of a split within the red-clump when investi¬ 
gating the luminosity function of the Bulge at latitudes |b| >5, along the minor 
axis. Soon enough. Me William and Zoccali ( 201 0[ > and |Nataf et al ( ]20 1 ()[ > provided 
a wider mapping of this split red-clump in the color magnitude diagram, provid¬ 
ing substantial evidence for the bright and faint red-clumps to be the consequence 
of having two over-densities of stars located at different distances, namely the two 
southern arms of an X-shaped structure both crossing the lines of sight. Detailed 
three-dimensional maps were later constructed by jSaito et al (2011]), based on the 
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use of red-clump stars observed in the near-IR survey 2MASS, which confirmed the 
suggestion of Me William and Zoccali | <j 201 Q|> that the Bulge is in fact X-shaped due 
to the prominent vertices of the B/P. Wegg and Gerhard (2013) modelled the distri¬ 
bution of red-clump stars, observed in Vista Variables in the Via Lactea (VVV) ESO 
public survey, providing the first complete mapping of the X-shaped Bulge. These 
X-shaped bulges are commonly observed in external edge-on galaxies and belong 
to the case of a pronounced B/P structure, that is simply the inner regions of the bar 
that grow out of the plane of the disc. 

Currently, the axial ratios of the bar are constrained to be about 1:0.4:0.3 with a 
bar size of about 3.1 —3.5 kpc diameter and the near end of the bar pointing towards 
positive Galactic longitudes. Until recently, the position angle of the bar was con¬ 
strained to a relatively large range of values between ~20—40 degrees with respect 
to the Suncentre line of sight. The uncertainty in the bar position angle is likely to be 
a consequence of measurements done across different latitudes in each study, thus 
finding a different position angle when looking at different distances from the Galac¬ 
tic plane. As a matter of fact, specific evidence for a longer flatter component of the 
bar, referred to as the Galactic long bar, has been presented in the literature based 


on near-IR star counts near the Galactic plane (Blitz and Spergel 1991, Stanek et al 

1994 

Dwek et al 1995 Binney et al 1997 

Bissantz and Gerhard 2002 Benjamin 

et al 

2005 Babusiaux and Gilmore 2005 

lattenbury et al 2007 Cao et al 

2013 

). 


This long bar is found to have an axis length of 4—4.5 kpc and ratios of 1:0.15:0.03. 
The position angle of this longer component has been constrained to ~45 degrees in 
such studies (e.g. Lopez-Corredoira et al 2007 Cabrera-Lavers et al 2007, Ham- 
mersley et al 2000 |Churchwell et al 2009 Amores et al 2013[ >. On the other 
hand, the recent model of the global distribution of red-clump stars from Weg g and] 
Gerhard ( |2013 i provided a precise measurement for the B/P Bulge position angle 
of 27 ± 2 degrees, in agreement with the studies done at larger distances from the 
Galactic plane. The nature of the long bar has been debated extensively in the lit¬ 
erature. Recently, Garzon and Lopez-Corredoira ( 2014[ > provided a theory where 
two co-existing bars, the long bar restricted to the plane latitudes and the B/P thick 
bar, could be present in the inner Galaxy. However, model observations of barred 


galaxies led [Martinez-Valpuesta and Gerhard (20111, Romero-Gomez et al (2011!, 
andjAthanassoula (2012) to strongly argue that the apparent long bar is an artefact 


associated with leading spiral features at the end of the shorter primary bar (the B/P 
Bulge). Furthermore, the co-existence of such independently large scale structures 
has not been seen in external galaxies. For this reason, the observed properties at¬ 
tributed to different bars in the Galaxy are more likely corresponding to a unique 
B/P bulge and bar structure formed by the buckling instability process. The long 
bar would then be explained by the interaction of the outer bar, with the adjacent 
spiral arm near the plane which produced leading ends that ultimately results in the 
measurement of a larger position angle (Martinez-Valpuesta and Gerhard, 2011| ). 

In the innermost regions (1 < 4, b < 2) the bar has been found to change its 
apparent inclination with respect to the line of sight which has been interpreted as 
evidence for a possible distinct smaller bar, referred to as a nuclear bar. However, 
models of a single bar (meaning those that do not include a distinct nuclear bar) 
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have also shown such a change in orientation in the inner regions, most likely due 
to the presence of a more axisymmetric concentration of stars in the central regions 
(Gerhard and Martinez-Valpuestaj 2012) >. 

Red-clump stars, although an excellent distance indicator for the bulge mean 
population, still suffer from the usual complications when looking towards the inner 
Galaxy such as disc contamination, extinction, and, to a minor extent, the effects 
of stellar populations. Furthermore, red-clump stars will map the distribution of the 
variety of stellar ages of the Bulge population, which is not defined beforehand. 
This uncertainties can be statistically handled when constructing and analysing the 
Bulge luminosity function, however their impact on the results will depend on the 
level of knowledge of the properties of each analysed field. For this reason. Variable 
stars, specifically RR-Lyrae, have recently provided a new perspective for the Bulge 
structural properties. Their well defined period-luminosity relation in the near-IR 
helps to overcome the effects of dust extinction in the inner Galaxy and they are well 
spatially distributed across the entire bulge. The period-luminosity relation makes 
RR-Lyrae an exquisitely accurate distance indicator that unequivocally traces the 
oldest Galactic population. 

Surprisingly, given the vast amount of different tracers that have confirmed a 
dominant barred structure, RR Lyrae have shown a remarkably different spatial dis¬ 
tribution compared to, for example, red-clump stars. While red-clump stars trace 
the position angle of the bar at all latitudes, a direct comparison with the RR Lyrae 
distance distribution provided by Dekany et al ( 2013| ) strongly suggests a different 
morphology for the oldest population in the inner Galaxy. Unlike the red-clump 
stars, the RR Lyrae stars show a more spheroidal, centrally concentrated distribu¬ 
tion. This structural component, populated by stars with ages larger than 10 Gyr, 
seems to be overlying with the B/P bulge. Figure [T] shows a comparison of the pro¬ 
jected mean distances obtained from RR Lyrae and those from the mean magnitude 
distribution of red-clump stars. The figure illustrates the structural difference be¬ 
tween the components traced by both distance tracers, with only red-clumps stars 
following the position angle of the bar. This result presented in Dekany et al ( 2013j > 
is perhaps the first purely morphological evidence suggesting a composite Bulge 
nature, with two different stellar populations overlapping in the inner Galaxy. 

The presence of more than one age/metallicity distribution within a B/P bulge has 
already been seen in dissipative collapse models (e.g. Samland and GerhardH2003j l 
and also in bulges from cosmological galaxy formation simulations (e.g. |Obreja| 
et al} 2013[). However, one must be very careful when further linking the different 


spatial distributions seen in the Galactic bulge with a distinct origin process, namely 
having a classical bulge and a secularly evolved B/P bulge that originated from 
the disc. Recently, [Ness et al (2014 1 gives caution to the fact that different spatial 
distributions and mean stellar ages can be found in pure B/P bulges without the need 
of a merger-origin structure to be present, as seen in an N-body + smoothed particle 
hydrodynamics simulation of a disc galaxy. Certainly, the fine details of the shape 
traced by RR-Lyrae will be achieved when the complete sample of RR-Lyrae from 
the VVV survey is available. The mapping of a wider area of the Bulge and the 
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Fig. 1 Upper panel: Spatial distribution in Galactic coordinates for the 7663 OGLE-III RRab stars 
in the bulge area of the VVV Survey, from [Dekany et al|j2013| l. The grey scale background shows 
the interstellar extinction map of |Schlegel et al[(1998) . The black rectangle denotes the region at 
b = —4 for which the mean distance of RR Lyrae and RC stars are compared in the lower panel. 
The lower panel shows the projected mean distances of RR Lyrae in black filled circles and of the 
red-clump stars as red open circles. Isodensity contours for the projected distance distribution of the 
RR Lyrae sample in the analysed latitude range are also shown. Mean distances of the red clump 
stars are from the mean magnitudes obtained in |Gonzalez et al| ( |2012| ) and calculated adopting 
an absolute magnitude of Mjc s ,rc = 1.71 mag. Distances for the sample of RR Lyrae have been 
presented in Dekany et al (2013 1 and were used here to derive the projected mean distance and Id 
width (black solid lines) to each line of sight. [Upper panel adapted from Dekany et al \2013\.\ 
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larger sample of sources available in each line of sight will allow for the calculation 
of precise mean distances with high spatial resolution. 


3 The age of the Milky Way bulge 


Right from the very early discovery of RR Lyrae towards the centre of the Galaxy 
( Baade| 1951} , perhaps the most common statement found in the literature address¬ 
ing the Galactic bulge stellar populations is: The Milky Way bulge stellar population 
is predominantly old. The reasons for such a statement are indeed very well funded 
and are described in this section. 

Without a doubt, deep photometric observations towards low extinction regions 
of the Bulge provided a defining view on the age of its stellar population thanks to 
the possibility of constructing colour-magnitude diagrams that reached the turn-off 
position - a useful indicator for the mean age of a given stellar population. However, 
the effects of differential reddening and, in particular, the uncertainty on the distance 
modulus of the stars towards a given line of sight does not allow the derivation of 
an absolute age estimation for the Bulge using the turn-off technique. To overcome 


these issues, Ortolani et al (19951 adopted a differential method, based on the po¬ 


sition of red-clump stars in the luminosity function. The very small dependence on 
age and metallicity for the mean magnitude of the red-clump allowed Ortolani et al 
( 1995} to match its position in the luminosity function to that of an old stellar popu¬ 
lation, the globular cluster NGC 6528, finding a remarkable agreement between the 
relative positions of the turn-off. This result provided strong support for the view of 
the Galactic Bulge being formed in a time scale shorter than 1 Gyr, thus making it 
as old as the globular cluster population. 

Dedicated photometric studies spread across other regions of the Bulge, and 
based on similar techniques further strengthen the conclusion of a stellar popula¬ 
tion with a mean age of ~ 10 Gyr, particularly setting a lower limit on ages higher 
than 5 Gyr (Zoccali et al 2003 1 Valenti et al 2013 Clarkson etal 201 1[ >. The major 
issue with these kind of studies is how to deal with the contamination of the fore¬ 
ground disc, in particular with the main sequence of the disc which lies right on top 
of the Bulge turn-off. Statistical decontamination methods, for example using disc 
control fields ( |Zoccali et al| [2003 ; Vale nti et al| [20 13), have been used to eliminate 
foreground stars to some extent. However, the contamination of foreground stars 
coupled with uncertainties in differential reddening and metallicity distribution ef¬ 
fects lead to uncertainties on age determination via the analysis of color-magnitude 
diagrams that remain of the order of 2 Gyr. 

Recently, Clarkson et al ( 2008} provided what is perhaps the cleanest, and thus 
most accurate, age determination of those studies based on the analysis of a colour- 
magnitude diagram (Fig. [2|. [Clarkson et al ( 2008f were able to decontaminate the 
turn-off position of the Bulge population from that of the nearby disc by using a 
mean proper motion criteria as the one shown by |Kuijken and Rich] ( j2002) > with 
a set of ACS WFC on HST observations taken over 123-orbit HST integrations 
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Fig. 2 Left panel: The metallicity of the microlensed dwarfs of the Bulge as a function of their ages 
taken from Bensby et al ( 2 013a}. Right pan el: The color magnitude diagram for proper motion- 
selected Bulge objects from Clarkson et al (2008), using similar mean proper motion criteria to 
Kuijken and Rich ( 2002) but with a 6(J detection requirement imposed. A set of isochrones with 
different metallicities and ages is overplotted to the color magnitude diagram. An alpha-enhanced, 
solar-metallicity isochrone at 11 Gyr represents the median sequence well above the turn-off. Also 
shown in the figure are sequences at metallicity [Fe/H]=(-1.009, -0.226, +0.491) and ages (8, 10, 
14) Gyr. [Left panel adapted from Fig. 15 in Bensby et al\\2013a\, ©ESO reproduced with per¬ 
mission. Right panel adapted from Fig. 20 in Clarkson et al\ \2008\ , ©AAS reproduced with per¬ 
mission. ] 


in the SWEEPS field (l,b)=(1.25, -2.65). The study of |Clarkson eTal) ( j20TT) later 
concluded that a fraction of only up to a 3.5% of the Bulge stars can be younger 
than 5 Gyr. Therefore, this study is in agreement with most of other studies where 
the Bulge population is found to be dominantly old. 

A large percentage of young stars is certainly not expected in a bulge that origi¬ 
nates purely from early dissipation or merging processes. These events would take 
place at early times in the evolution of the galaxy and would occur rapidly. Al¬ 
though some intermediate-age stars could later be added to the bulge by diffusion 
from the inner disc, the bulk of the bulge stellar population would be old. On the 
other hand, a significant population of young stars could be found — and might be 
actually expected — in bulges formed via disc instabilities. So if the Milky Way 
bulge is certainly dominated by a component formed via disc instability, as seen by 
its structural properties, then where are those young stars? It is only very recently 
that a possible answer to this question has been brought to the table, and it was 
using a completely different approach: the high-resolution spectroscopy of Bulge 
microlensed stars. The microlensing event, which produces the brightness magni¬ 
fication of a dwarf star in the bulge due to the passing of a foreground lens star, 
provides a unique opportunity to obtain high resolution spectra of this otherwise un¬ 
reachable target. Bensby et al ( 2013a] ) have collected enough micro-lensed dwarfs 
to investigate the overall metallicity distribution and also their age. While the metal- 
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licity distribution of the microlensed dwarfs has been found to be in good agreement 
with that of the bulge giant stars, their age distribution presented a significant num¬ 
ber of young stars: nearly 22% of the micro-lensed dwarfs were found to be younger 
than 5 Gyr. As the number of analysed micro-lensed stars increases, these findings 
will be further confirmed or disproved with better statistics in order to refine the 
actual percentages of young and old stars found in the bulge. Particularly, notice the 
disagreement between the 22% of stars younger than 5 Gyr found in the microlens- 
ing sample and the corresponding 3.5% fraction of young stars given in the work of 
|Clarkson et al| ( f20TT) >. 


4 The chemical abundances of the Milky Way bulge 
4.1 The metallicity distribution 


The characterisation of the chemical abundance of Bulge stars is perhaps the field 
that has evolved most quickly thanks to the advent of multi-object spectrographs in 
large telescopes. Different surveys have — and still are — pointing towards different 
regions of the bulge and collecting samples of thousands of stars. The argument for 
such surveys is clear: it became evident that a few low extinction regions were no 
longer representative of the global chemical abundance patterns of the Bulge. First 
attempts to derive the metallicity distribution of the Bulge based on low resolution 
spectra (e.g. Sadler et al 1996[ Ramirez et al 2000[>, together with a sm all number 
of available spectra obtained with high resolution (McWilliam and Rich |1994| Ful-| 
bright et al 2006), had shown from the start that the mean population of the bulge 


was overall metal-rich. The shape of the metallicity distribution was, on the other 
hand, less clear. Zoccali et al ( 2003| ) used a set of observations using WFI to obtain 
a photometric metallicity distribution based on the colour of red giant stars. They 
found a rather broad metallicity distribution, in good agreement with that derived 
from spectroscopy in the same field but with a larger statistical sample, with [Fe/H] 
values ranging from -1.0 to 0.4 and which peak at solar metallicity. 

With a well-characterised metallicity distribution in Baade’s window, it was time 


to answer the following question: how spatially uniform were these properties? Min- 


niti et al ( 1995| had already discussed the possibility for a metallicity gradient in the 


Bulge, impressively enough based on low resolution spectra of less than a hundred 
giant stars. It then became the era of multi-object spectroscopy where hundreds of 
stars could be observed in one single shot. Using FLAMES on the Very Large Tele¬ 
scope, Zoccali et al (2008) derived the metallicity distribution for different fields 
along the Bulge minor axis at different latitudes (b = —4, —6, and —12). They 
found a clear metallicity gradient of ~0.6 dex/kpc, with mean metallicities vary¬ 
ing from -0.4 at the largest latitudes and up to solar metallicity at b=-4. This gradi¬ 
ent has since then been confirmed thanks to several subsequent observations across 
different regions and the variation in the metallicity distributions have been further 



































10 


Oscar A. Gonzalez and Dimitri A. Gadotti 


characterised (Johnson et alj 2011 

Uttenthaler et alj 2012|[Johnson et al[|2013||Ness| 

et al 2013aJ Rojas-Arriagada et a 

2014). 



1 
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Galactic Longitude (1) 


Fig. 3 Left panel: Metallicity distributions for a compilation of studies taken from 


Johnson et al 


\2() 13} . The (1,b) = (+1,-4), (0, —6) and (0,-12) fields are fro m|Zoccali et al | j2008), th e (0,-8) 
field is from |Johnson et al] \20 1 I } , the (+5,-3) fie ld is from | Gonzalez ct al| ( |201 1 } , and the 
(—5.5, —7), (—4, —9), and (+8.5,+9) fields are from |Johnson et al||20l 3} . Right panel: Map of 
the mean photometric metallicities of the bulge constructed with the VVV survey data from |Gon-| 
|zalez et al|j2013} . [Left panel adapted from Fig. 8 in Johnson et aI [2013), ©AAS reproduced with 
permission. Right panel adapted from Fig. 2 in \Gonzalez et al\\2013\ , ©ESO reproduced with 
permission .] 


Gonzalez et al (|2013|l recently complemented these results by presenting a pho¬ 


tometric metallicity map, constructed with the same technique used by Zoccali et al 
( 2008[ > but based on the Vista Variables in the Via Lactea (VVV) ESO public sur¬ 
vey, for almost the entire Bulge region providing the global picture of the Bulge 
metallicity gradient. The metallicity gradient is therefore strongly established by 
an increasing number of spectroscopic studies obtained with different techniques 
and stellar samples (see Fig|3]for a compilation of the latest results). However, the 
metallicity distributions obtained in the innermost regions of the Bulge (\h\ < 4) 
based solely on high-resolution, near-infrared spectroscopy have provided evidence 
for the flattening-out of the gradient in the inner 700 pc (jRarmrez et al| |2000| Rich 
[etaIl|20T2] >. 


What is the implication of finding such a metallicity gradient in the Bulge? At 
first, similarly to the domination of old ages found in Bulge stars, the metallic¬ 
ity gradient was interpreted as direct evidence for a bulge formed as a classical 
bulge via mergers in the early stages of the galaxy, similarly to elliptical galaxies. 
It was also interpreted as evidence against the secular evolution scenario, since it 
was thought that bars would mix the stellar orbits well enough to erase any existing 
vertical gradient. Models of bar formation in disc galaxies, however, proved oth¬ 
erwise, showing that a bar might produce a gradient similar to the one seen in the 
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Milky Way depending on the original disc radial gradients (|Martinez-Valpuesta and 


Gerhard 2013), vertical gradients (Bekki and Tsujimoto 2011), or both (Di Matteo 
et al||2014). ' 


However, the existence of metallicity gradients has also been interpreted as a 
consequence of having two or more underlying components each one with a char¬ 
acteristic metallicity distribution. This mixing of components would naturally pro¬ 
duce a variation on the mean metallicity according to the bulge region which is 
being studied. Evidence for such a multiple component scenario has been suggested 
based on a bimodal metallicity distribution of red clump stars in Baade’s window 
by Hill et al ( 2011} and has been also suggested from a similar bi-modality seen in 
the metallicity distribution of microlensed dwarfs. Both of this distributions show 
a metal-poor and a metal-rich peaks located approximately at [Fe/H] ~ —0.3 and 
[Fe/H] ~ +0.3, respectively. Recently, the same bi-modality has also been found 
in the metallicity distributions based on the Gaia-ESO survey observations of the 
Bulge (Rojas-Arriagada et al 2014) . By producing a Gaussian decomposition of 
the metallicity distribution functions, Rojas-Arriagada et al (2014) showed a clear 
bi-modality in all the analysed fields with relative sizes of components depending of 
the specific position on the sky. This change in the relative sizes of each component 
can be clearly seen in Fig. [4] 

On the other hand, the ARGOS Galactic bulge survey (Freeman et al 2013), 
which consists of the largest sample of homogeneously analysed RC stars, con¬ 
structed large-number statistics metallicity distributions at different latitude stripes, 
using a total of more than 10,500 stars located within a galacto-centric radius of 3.5 
kpc. The overall metallicity distribution of the ARGOS survey was interpreted in 
Ness et al (2013a)), as being composed of five Gaussian components. Each of these 


components would be sampling a different stellar population and thus any changes 
in their relative contribution fraction as a function of latitude could be the origin 
of the observed mean metallicity gradient seen in the Bulge. The metallicity distri¬ 
butions for stars within Galactic longitudes l = ±15 and latitudes b = —5, —7.5, 
and —10 from the ARGOS survey Ness et al ( 2013a) are also shown in Fig. [4] This 
effect has lead the ARGOS survey to suggest the three main components of the 
metallicity distribution to be associated with the metal rich B/P bulge (mean [Fe/H] 
~ +0.15), the thick B/P bulge (mean [Fe/H]~ —0.25) and the inner thick disc (mean 
[Fe/H]~ -0.70). 

Although the results from Ness et al (2013a i are based on the identification of 
several components, it can be safely understood that most studies converge into a 
similar conclusion, the bulge metallicity distribution is the result of not a single 
but of a mixture of populations with at least two main components, one metal-poor 
and the other metal-rich. Certainly, this is everything that can be concluded from 
the metallicities alone and one has to be extremely cautious when attempting to 
link these different components with a given bulge formation scenario. All evidence 
needs to be considered when interpreting results in terms of such scenarios, i.e. such 
an assessment must also include the kinematics, spatial distribution and ages of the 
stars. 
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[Fe/H] [Fe/H] [Fe/H] 


Fig. 4 Upper panels: Metallicity distributions for the five bulge fields from the ESO Gaia survey, 
presented in |Rojas-Arriagada et al|f2014f . Black dashed and solid lines show the components iden¬ 
tified in each field as dashed Gaussian functions, with the sum of them shown as a solid line. The 
three minor-axis fields located at (+1,-4), (0,-6), and (-1,-10) are shown in the central panels, and 
the lateral fields (+7,-9) and (-10,-8) at the left and right. Lower panels: Metallicity distributions 
from the ARGOS survey from |Ness et al| ( {2013b[ >. From left to right at b = —5, —7.5 and —10, for 
1 = ±15. The different contribution of the adopted Gaussian components are marked in each field, 
with the three main components being A, B and C. [Upper panel adapted from Fig. 6 in \Rojas^\ 
Arriagada et al {2014) . ©ESO reproduced with permission. Bottom panel adapted from Fig. 1 in 
Nessjtt alt201 Mu | 


4.2 The Bulge alpha-element abundances 


The a-element abundances can provide us with further constraints for the origin of 
the Bulge stellar populations, specifically with respect to its formation time-scale. 


Tinsley (1979) suggested that the ratio of [a/Fe] compared to [Fe/H] is a function of 


the time delay between SNe II, which produce both a- and iron-peak elements (e.g. 
[Woosley and Weaver||1995|l, and SNe la , which yield mostly iron-peak with little a- 
element production (e.g. Nomoto et al 1984). Therefore, only after sufficient time 


has passed for the SNe la events occur, the [a/Fe] ratio will decline from the SNe 
II value. Clearly, the critical ingredient on this relation is the SNe la delay time, for 
which different production channels might be present. 

In the Bulge, the a-element abundances of Bulge stars with [Fe/H] <-0.3 have 
been found to be enhanced over iron by [a/Fe]~+0.3 dex (McWilliam and Rich 
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1994 Rich and Origlia 

2005, 

Cunha and Smith) 2006; |Fulbright et al| 2007J 

Lecureur et al 2007 Rich et al 

2007) calling for a fast formation scenario, while 


metal-rich stars [Fe/H]>-0.3 showed a decrease in [a/Fe] reaching solar values for 
metallicities larger than Solar. However, is important to note that not all elements 
were found to follow the same yield trends. 

A relative approach has been commonly adopted instead of an absolute inter¬ 
pretation of the Bulge a-element ratio. The direct comparison of [a/Fe] values in 
Bulge stars against those of other galactic components then provides a relative time 
constraint on the bulge formation. Fulbright et al ( 2007) ), Zoccali et al ( 2006] ), and 
|Lecureur et al| ( |2007] ) all came to the conclusion that the [a/Fe] ratio was enhanced 
by nearly +0.1 dex with respect to the trends of both the thin and the thick disc, thus 
implying a shorter formation time scale for the bulge than from both discs. However, 


as first pointed out by Melendez et al (2008), the bulge a-element over-enhancement 
with respect to the thick disc was a result of systematic offsets between abundance 
measurements in dwarf stars of the disc and giant stars from the Bulge. Later on, 
Alves-Brito et al (2010) and Gonzalez et al (201 1]) confirmed that when giants from 


the both bulge and the disc are homogeneously analysed the Bulge followed the 
same over-abundance in a-elements as the thick disc, both being enhanced with re¬ 
spect to the thin disc at metallicities [Fe/H]<-0.2. At solar metallicities, the Bulge 
stars are found to be a-poor, as poor as those of the thin disc. The way these trends 
are interpreted is that the metal-poor population of the bulge underwent a similarly 
fast formation scenario to that of the thick disc, while the metal-rich population of 
the bulge must have had a longer formation time scale, in similar time-scale to that 
of the thin disc stars. Similar conclusions have been reached in several other studies 
carried in different regions of the Bulge (Bensby et al 2010 201 1[ Ryde et al 2010 


|Hill et al| [20TT| |Johnson eTd||20TTi|2013||2014). 

However, open questions remain regarding the chemical similarities of Bulge 
stars with those of the thick disc stars, particularly in light of a few recent findings. 


Bensby et al (2013a I suggested that the position in the [a/Fe] - [Fe/H] plot where 


[a/Fe] starts to decrease (referred to as the knee in the literature) is located at higher 
metallicities in the Bulge than in the thick disc. The position of the knee in the bulge 
may be 0.1-0.2 dex higher in metallicity in the Bulge than in the thick disc thus 
suggesting that the chemical enrichment of the metal-poor bulge has been some¬ 
what faster than what is observed for the local thick disk. As the sample of Bulge 
micro-lensed dwarfs increases, it would be of great interest to further confirm the 
findings of Bensby et al ( 2013a| ). As a matter of fact, a similar result was proposed 
by Johnson et al ( 2014] ) who also added the analysis of Fe-peak elements finding in 
particular that Co, Ni, and Cu appear enhanced compared to the disc. It is impor¬ 
tant to recall that the results presented in Johnson et al ( |2014| > have been obtained 
by comparing Bulge giants to dwarf stars from the local disc. This technique has 
been shown to suffer from systematic offsets Melendez et al ( 2008] ). However, the 
detailed analysis by Johnson et al ( 2014| ) has been carefully calibrated internally so 
it would be of great interest to confirm if these results are also found when bulge gi¬ 
ant stars are compared to (inner) disk giant stars. These findings certainly highlight 
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the importance of the future multi-object spectroscopic surveys on different galactic 
components to obtain a definitive answer. 


5 The Bulge Kinematics 


The observational properties of the Bulge regarding the morphology, age and chem¬ 
ical abundances seem to be independently providing evidence for a rather complex 
bulge stellar population, where at least a metal-poor, alpha-enhanced population of 
old stars co-exists with a metal-rich, alpha-poor population of both old and a frac¬ 
tion of young stars (~ 22% according to Bensby et al ( |2013a| > but see|Clarkson et al| 
poll) ). ] t is then natural to evaluate how the kinematics of each of these popula¬ 
tions can help to solve the puzzle and answer the question on what is the nature of 
each of these components. 


Following the initial attempts to constrain the kinematics of the bulge stars by 


measuring their radial velocities (Frogel and Whitford 1987 Rich 1988 

1990; 

Terndrup et al 1995( Minniti 1996} Sadler et al 

1996 Tiede and Terndrup 

1997) 


our understanding of the dynamical characteristics of the Bulge has gained an out¬ 
standing level of advancement thanks to the recent spectroscopic surveys that are 
able to sample radial velocities of thousands of M giants and red-clump stars across 
the Bulge. 

In the most general view of Galactic bulge kinematics, the bulge is known to lie 
between a purely rotating system and a hotter system supported by velocity disper¬ 
sion, with a V max /a = 0.65 (Minniti and Zoccali 2008 . The Bulge Radial Velocity 
Assay survey (BRAVA; Howard et al 2008[ Kunder et al 2012| > was the first to 
provide a broad view of the bulge kinematics, perhaps finally allowing us to start 
looking at the Bulge from an extragalactic perspective. The BRAVA survey pre¬ 
sented the mean radial velocity and velocity dispersion as a function of longitude, 
at different Bulge latitudes, showing evidence for cylindrical rotation of the bulge 
(BRAVA; |Howard et al[|2008}|Kunder et al||2012[ > which is a characteristic feature 
of box/peanuts originating from secularly evolved bars. 

The cylindrical rotation seen in M giants of BRAVA, later confirmed using red- 
clump stars by the ARGOS (Ness et al 2013b| > and GIBS ( Zoccali et al| 2014| > 
surveys, was modelled by Shen et al| ( |2010| ), who conclude that the bulge could 
not have any more than 8% of the disc mass in the form of a classical spheroid to 
reproduce BRAVA observations. Figure[5]shows the agreement between the BRAVA 
survey measurements at different latitudes compared to the pure disc models from 
Shen et al (2010 1 . Also shown in Fig.[5]are the radial velocity and velocity dispersion 


maps for the Bulge constructed by Zoccali et al ( 2014| based on the GIBS survey 
observations. 

However, dynamical models such as the one presented by Saha et al ( j2012| > have 
shown that if a spheroidal component, i.e. a classical bulge, was already present 
when the bar was formed then the classical bulge could spin-up and rotate faster 
than expected for its dispersion supported nature due to the effects of the bar poten- 
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Galactic Longitude (degrees) 


Galactic longitude (deg) 


Fig. 5 Left panels: Velocity dispersion profile (bottom) and rotation curve (top) for latitudes b = 
-4, -6, and -8 from Kunder et al (2012 1 . Solid lines represent the models from Shen et al 1 2010) . 
Right panels: Radial velocity (top) and radial velocity dispersion (bottom) surface in the longitude- 
latitude plane constructed from the measured rotation profiles at negative latitudes from the GIBS 
survey by |Zoccali et al|{2014[ l. Grey points show the positions of the observed fields by the survey, 
while the black contour lines are labeled with the relevant velocity dispersion in km/s. [ Left panel 
adapted from Fig. 11 in \Kunder et al\{2012 , ©TAS reproduced with permission. Right panels 
adapted from Fig. 10 and 11 in\Zoccali et a 1 1 2014} , ©ESO reproduced with permission.] 


tial ( jSaha and Gerhard 20131. In these conditions, the detection of such a compo¬ 
nent based in kinematics alone would be very difficult ( |Gardner et al||2014] l. Indeed, 
in order to understand the nature of these components the analysis of a connection 
between the kinematics and other stellar properties such as metallicity seems to be 
a key factor. Babusiaux et al ( 2010| investigated the connections of metallicity and 
kinematics, the latter based on radial velocities and proper-motions, for the sample 
of |Zoccali et al| ( |2()08| > at different latitudes along the bulge minor axis. They found 
that the high metallicity stars ([Fe/H]>-0.25) show a larger vertex deviations of the 
velocity ellipsoid than their metal-poor ([Fe/H]<-0.25) counterpart. Furthermore, 
metal-rich stars showed an increase in their velocity dispersion with decreasing lat¬ 
itude (moving closer to the galactic plane), while metal-poor stars show no changes 
in the velocity dispersion profiles. This information led Babusiaux et al ( 2Q10j > to as¬ 
sociate the more metal-rich stars with a barred population and the metal-poor stars 
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with a spheroidal component or even the inner thick disc. The rotation curves and 
dispersion profiles of the large sample of stars from the ARGUS survey (Fig. [6} led 
|Ness et aT](|2013b]l to reach a similar conclusion to that of|Babusiaux et al|(|2Q10|>. 



Fig. 6 Rotation (top panel) and velocity dispersion (bottom panel) across the bulge for the 16600 
stars from the ARGO survey with [Fe/H] > -1.0 within < 3.5 kpc of the Galactic Centre from 
|Ness et a?|(2013b\ The three plots correspond to different metallicity bins, from left to right in 
decreasing [Fe/H]. Note that the discrete bins are used to represent stars of components A, B and 
C from left to right shown in Fig[4] Although the rotation curves are similar, the dispersion clearly 
demonstrates the difference in kinematics of stars with [Fe/H]> -0.5 and with [Fe/H]< -0.5. There 
are 3100, 8600 and 4900 stars in each plot, from left to right. The red diamonds are b = -5, the 
yellow stars are b = -7.5, the blue rectangles are b = -10 and the black circles are b = +10. [Figure 
adapted from Fig. 6 in\Ness et al\{2013tty.] 


An additional piece of the puzzle has been provided thanks to the recent develop¬ 
ment of our understanding of the structural properties of the Bulge, in particular the 
discovery of the X-shape. Ness et al (2013b i and Vasquez et al ( 2013| > investigated 
in more detail the connection between the X-shape bulge and its chemo-dynamical 
properties. Both studies showed that only the metal-rich stars ([Fe/H]> —0.5) trace 
the split red-clump in the luminosity function and therefore belong to the X-shaped 
bulge, while the metal-poor stars do not share the same split in magnitude. 

Finally, important constraints on the bulge formation history can be obtained by 
looking at the radial velocities of a large number of stars. However, also having 
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Fig. 7 Comparison of the kinematics of the bright and faint red clump of the Bulge with respect 
to the kinematical model for a strong boxy-peanut bulge {Debattista et al 2005) obtained from 
|Vasquez et al (2013 1 . From the model, two samples were selected from the two overdensities 
formed by the near (red) and far (blue) arms of the boxy-peanut stellar distribution in the line of 
sight for (l,b) = (0, -6). Colour lines over U, V, and W histograms correspond to the median value 
for each distribution. [Figure adapted from Fig. 7 in \Vdsquez et al\$201 jp , ©ESO reproduced with 
permission.] 
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information on the proper motions for some of those stars in order to reconstruct 
the three-dimensional characteristics of the stellar orbits may become fundamen¬ 
tal to disentangle the more complex characteristics of the Bulge. Although proper 
motions can be derived for large samples of stars based on photometric informa¬ 
tion obtained with relatively short exposure times, they require a long time baseline 
(~ 10 yr). Therefore, to obtain a very precise astrometric solution in order to reach 
the required accuracy of a few mas/yr across the bulge is perhaps more a luxury than 
a requirement. Indeed, studies that have used proper motion information in specific 
fields, combined with radial velocity measurements, have been able to unravel the 
perhaps otherwise hidden complexity of the Bulge (e.g. Zhao et al 1994[ Soto et al 


2007 1 Babusiaux et al| 2010[>. In this context, Vasquez et al (2013 provided for the 


first time an analysis based on both the radial velocities and proper motions for both 
arms of the X-shape bulge, and were thus able to derive the complete space veloci- 
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ties in the U, V, W Galactic Cartesian system for a sample of spectroscopic targets 
in the field (l,b)=(0,-6). As shown in Fig.[7] |Vasquez et~al ( |2013 i showed how the 
closer over-density of the X-shaped bulge shows an excess of stars moving towards 
the Sun, and the far overdensity shows an excess of stars receding from the Sun, as 
expected from stars on elongated orbits streaming along the arms of the X-shaped 
bulge. A wider mapping of these are the key signatures of the detailed kinematics 
of the bulge, thus expanding the study of Vasquez et al (2013 ) to other Bulge re¬ 
gions, will become an important ingredient for an accurate modelling of the bulge 
and will perhaps allow a better identification of its different components [Gardner] 
[etaH ( |201^] ). 


6 The Milky Way bulge in a nutshell 

The properties of the Galactic bulge described above can be summarised as follows: 

• The Bulge stellar population properties show independent evidence for a multiple 
component scenario, with different morphological and dynamical characteristics. 
The metallicity distribution can be separated in, at least, a main metal-poor and a 
metal-rich components. Their different contributions across the Bulge produces 
a vertical metallicity gradient of ~0.6 dex/kpc which seems to flatten in the inner 
700 pc. From a morphological/structural viewpoint, there is strong evidence of 
a box/peanut, while more recent work points out the existence of a component 
with different geometrical properties, as seen in the spatial distribution of RR 
Lyrae stars, which trace the oldest Galactic population (Fig. [T}. The next two 
bullet points connect these structural components to differences observed also in 
chemical content, age and dynamics. 

• There is a metal-rich (mean [Fe/H ]~0.3), alpha-poor population of stars in the 
bulge, composed mostly by old stars but with a fraction of young stars. This 
population of stars rotates cylindrically and shows a large vertex deviation con¬ 
sistent with the bar structure traced with a position angle of ^27 deg. At latitudes 
\b\ > 5 the inner parts of the bar have grown out of the disc plane, originating 
the box/peanut. These metal-rich stars further show the split red-clump in the 
luminosity function tracing the X-shape of the Bulge. 

• The metal-poor bulge population (mean [Fe/H |—-0.3) is composed predomi¬ 
nantly by old stars and it shows an alpha-enhancement similar to that of the local 
thick disc. The kinematics of these stars follow a more spheroidally distributed 
population than the one traced by the metal-rich stars, consistent with the struc¬ 
ture traced by the bulge RR Lyrae stars. Furthermore, these stars do not trace the 
X-shape morphology of the Bulge. 


















The Milky Way Bulge: Observed properties and a comparison to external galaxies 


19 


7 The Milky Way bulge in the context of external galaxies 
7.1 The X-shaped Bulge of the Milky Way: how rare is this? 


It can be argued that while in the context of Galactic research box/peanuts are 
somewhat a recent discussion, in Extragalactic studies such structures are known 
for about twice as much the time. So it came as no surprise to the Extragalactic 
community when evidence suggested that the Milky Way has a box/peanut, espe¬ 
cially because there is also evidence that it has a bar. Perhaps the first mention in 
the literature about these deceptively unusually-looking structures is from Burbidge] 
( 1959| >, referring to the prototypical example that is NGC 128 (see 
1961[>. More detailed investigation came with|de Vaucouleurs (19741, 


and Burbidge 


also 


Sandage 


Jarvis ( 1986|> and 
Dos Anjos fl 987 


Shaw (1987). But it was in the pioneering studies of de Souza and 


that the major step of connecting these structures with bars was 
made for the first time from an observational viewpoint, using a statistical argument. 
Basically, they argued, the frequency of box/peanuts in edge-on lenticulars is simi¬ 
lar to that of bars in face-on lenticulars, consistent with the idea that box/peanuts are 
bars seen at a different projection. This conclusion was corroborated years later by 
Liitticke et al (2000), who reported a fraction > 40 per cent of box/peanuts in disc 


galaxies covering most of the Hubble sequence (from SO to Sd classes). 

The complicating factor here is of course the fact that bars are difficult to be 
seen when the inclination angle of the galaxy is too large. Therefore, simulations of 
barred galaxies played a major role here. In fact, the starting point for this observa¬ 


tional connection between bars and box/peanuts was the work published in Combes 


and Sanders|( jl98l) . These authors have shown, using collisionless simulations, that 
bars seen at a given edge-on projection show a very characteristic peanut-like mor¬ 
phology. 

A number of studies came thereafter dedicated to extend this connection into a 
dynamical context. Kuijken and Merrifield ( 1995) came up with an ingenious di¬ 
agnostic to test in this context whether box/peanuts are just bars seen at a different 
projection. This consisted in producing diagrams in which the line of sight velocity 
is plotted against the galactocentric radius for highly inclined or edge-on systems. 
By producing such diagrams corresponding to orbits in a purely axisymmetric po¬ 
tential, and orbits in a barred potential, they showed that the presence of a bar pro¬ 
duces a clear distinctive signature. Because at a region around the bar corotation 
radius (the radius at which the pattern speed of the bar matches the local circular 
speed) there are no close, non-self-intersecting orbits available, clear gaps appear 
in this diagram, producing a figure-of-eight pattern. The matter became then just 
to produce such diagrams for galaxies presenting box/peanuts in order to test for 
the presence of a bar. Kuijken and Merrifield ( 1995) did that for NGC 5746 and 
NGC 5965, providing observational evidence that box/peanuts and bars are related 
phenomena. 

More evidence was produced in Merrifield and Kuijken (|1999 1 and Bureau and 


Freeman (1999). To this point, almost all galaxies with box/peanuts studied showed 
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evidence of a bar. In only a few extreme cases the box/peanut could have formed 
through accretion of external material. In addition, none of the galaxies without 
box/peanuts showed signatures of a bar. Further development also happened in the 
theoretical background. Bureau and Athanassoula (|1999[) refined and corroborated 


ten and Merrifield (1995), while Athanassoula and Bureau 


the orbital study of Kuij 
(1999) provided strong support to the bar detection diagnostic with hydrodynamical 
simulations. 

It must be noted that the detailed morphology of box/peanuts, i.e. if they either 
have a boxy shape, a peanut shape or an X shape, depends on projection effects, 
as well as the strength of the box/peanut. More about X shape bulges can be found 
in Laurikainen & Salo (this volume), but essentially the X shape is more clear in 
the strongest peanuts. Techniques such as unsharp masking are able to reveal X 
shapes more clearly when the peanut is not so pronounced. Therefore, all these 
morphologies come actually from the same physical process, i.e. the increase in the 
vertical extent of stellar orbits in the inner parts of the bar. Both simulations and 
observations point out that box/peanuts extend to galactocentric distances which 
are about a third to a half of the bar semi-major axis (see Erwin and Debattistaj 
2013). What causes this change in the orbits vertical extent is reviewed in detail by 


Athanassoula (this volume). In the Milky Way, we see the extension of the B/P to 
galactocentric distances of ~ 1.5 kpc, which is nearly two-thirds the length of the 
semi-major axis of the long bar of 2-2.3 kpc (Hammersley et al 2000 1 Benjamin 
et al 2005 1 Lopez-Corredoira et al 2007| l. 


A review on the observed properties of box/peanuts in external galaxies is given 
by Laurikainen & Salo (this volume). They also discuss a structure called barlens, 
which is interpreted as the projection of box/peanuts when seen face-on. This struc¬ 
ture was noticed by Laurikainen et al (20051, who included a model to fit barlenses 
in their image decompositions. Later, based on Fourier analysis, Laurikainen et al 
( j2007f suggested that barlenses are part of the bar, while |Gadotti| ( j200'8| also noticed 
their existence in a sample of local barred galaxies. However, only in Laurikainen] 


et al (20111 the term ’barlens’ is introduced as a new morphological feature in galax¬ 


ies. Very recently, more detailed studies have made a robust connection between 
barlenses and box/peanuts (Laurikaine n et al]|2014[ [Athanassoula et al 20141. This 
connection implies that also on the plane of the disc, the stellar orbits in the in¬ 
ner part of the bar become wider. Ligure [8] describes schematically the connection 
between bars, box/peanuts and barlenses. 

The presence of barlenses in external galaxies suggests of course that our own 
Milky Way may have such structural component. Since the stars in barlenses seem 
to be contained within the disc plane, in the Milky Way, they are seen in projection, 
in the foreground and background with respect to the Galactic centre. They thus 
complicate even further the interpretation of observations of the Milky Way, such as 
those discussed above. These stars are stars within the bar and then have chemical 
properties and ages similar to those of other bar stars - they are box/peanut stars. 
Lrom a kinematical point of view, barlenses are different from discs and this is a 
promising avenue to separate them from the other stellar populations seen from the 
Sun at the direction of the Milky Way central regions. 
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Fig. 8 Connection between bars, box/peanuts and barlenses. The diagram on the left shows a 
schematic representation of a barred galaxy seen a) face-on, and b) edge-on. Box/peanuts can 
be seen at edge-on projections, but the flat, more extended part of bar is difficult to realize, as 
its vertical extent is similar to that of the disc. The same galaxy seen face-on would reveal the 
bar and barlens. The barlens and the box/peanut appear to be the same structure seen at different 
projections [see Laurikainen & Salo (this volume) and Athanassoula (this volume)]. The panels on 
the right show an R-band image of NGC 4608 (top), where the bar and barlens can be clearly seen 
(see also Laurik ainen et al| |2005 2007 2011). The bottom panel shows a residual image, after 
the subtraction of a 2D model of the bulge, bar and disc of this galaxy. In the residual image the 
barlens stands out even more clearly. The red circle points out the barlens. This circle was not in 
the original Gadotti (2008) paper, but added now that we understand that the structure is a barlens. 
The red arrows point out empty regions in the disc within the bar radius, where stars from the disc 
were captured by the bar. [Right panels adapted from\Gadotti\^2008\.] 


7.2 Bulge formation scenarios 


The formation of bulges in general, and of the Milky Way bulge in particular, have 


been discussed many times elsewhere (e.g. Aguerri et al 2001, Kormendy and Ken- 

nicutt 2004; Athanassoula 

2005 

Laurikainen et al 2007 Hopkins et al 2010 

Fisher and Drory 2010 Gadotti 

2012 see also Bournaud, and Brooks & Chris- 


tensen, this volume). Here we will assess each bulge formation scenario in the light 
of evidence obtained from data on the Milky Way bulge, as presented above. There 
is mounting evidence that the Milky Way has a bar and a box/peanut, and thus 
secular evolution processes induced by bars in disc galaxies must have played a 
non-negligible role in the evolution of the Galaxy. On the other hand, a scenario in 
which the Galactic bulge was formed in violent processes such as mergers has weak 
support from data. In fact, an important question that observers must focus now is 
whether the Galaxy has a merger-built central stellar component at all. 





















































22 


Oscar A. Gonzalez and Dimitri A. Gadotti 


7.2.1 Bulges formed via disc instabilities 


Dynamical disc instabilities can originate bars, spiral arms and ovals, the latter being 
just a distortion in the disc stellar orbits that make them acquired less circular orbits, 
but not as eccentric as those in bars. All these structures, being non-axisymmetric, 
produce perturbations in the galaxy potential, with the result that material (gas and 
stars) within the corotation resonance radius loses angular momentum, whereas ma¬ 
terial outside the corotation radius absorb this angular momentum. The effect is par¬ 
ticularly important for the collisional gas component, which thus falls towards the 
centre. At some point, the in-falling gas gets compressed and form stars, contribut¬ 
ing to a rejuvenation of the stellar population in the central regions (e.g. |Gadotti and| 
|dos Anjosl|2001||Fisher||2006||Coelho and Gadotti[|201 l||Ellison et al||201 l|i. 

Because most disc galaxies with bars should have one or two Inner Lindblad 
Resonances (ILRs) near the centre (~ a hundred to a few hundred parsecs from it), 
the in-falling gas cannot reach the galaxy centre immediately (see Fig.|9|. Instead, it 
usually forms an inner disc, decoupled from the large-scale disc. These inner discs 
may form nuclear bars and spiral arms, as often observed, and are often called disc¬ 
like bulges (see e.g. Athanassoula 2005 1 Gadotti 2012} or discy pseudobulges, to 
contrast with the fact that box/peanuts are often called as well pseudobulges (see 
Kormendy and Kennicutt 2004} . At the ILRs gas can get accumulated and com¬ 


pressed, often forming a star-bursting nuclear ring. 



Fig. 9 Angular speed of stars in circular orbits in a potential reproducing a disc galaxy with a bar, 
as a function of the galactocentric distance (in arbitrary units). The bar pattern speed is represented 
by the solid horizontal line (Qb)■ The epicyclic frequency of the stellar orbits is denoted by K. 
Barred galaxies present several dynamical resonances. This figure shows schematically how four 
of the main resonances come to be. Whenever the bar pattern speed is equal to Q or Q ± k/2 this 
is the position of a dynamical resonance. From the centre outwards, the resonances depicted here 
are: the Inner Inner Lindblad Resonance, the Outer Inner Lindblad Resonance, Corotation, and the 
Outer Lindblad Resonance. In this case, the main families of stellar orbits change their orientation 
by 90 degrees at each resonance. This effect is at the origin of a number of dynamical effects in 
barred galaxies, in particular the transfer of angular momentum from material inside the corotation 
radius to material outside this radius, and the resulting formation of disc-like bulges. 
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However, it must be noted that simulations indicate that, initially, the infall of 
gas within a bar corotation radius occurs rapidly after the formation of the bar, with 
a time scale of the ord er of 10 8 years, i.e. a dynamical time ( Athanassoula[ 1992 


Emsellem et al 


2014)0 This means that disc-like bulges not necessarily present 


ongoing star formation or a very young stellar population, if the bar has formed 
long ago and have been able to push most of the gas to the centre quickly, and if 
the gas content in the disc is not being replenished. Sheth et al|( 2Q08}|2012] i present 
results that suggest that the first long-standing bars ^formed after redshift 1. These 
bars are still the ones seen at redshift zero, since bars are difficult to destroy, unless 
the disc is extremely gas-rich (see|Athanassoula et al 2005 Bournaud and Combes 


2002 [ [Bournaud et~aI||2005)|Kraljic et al| 2012| i. This means that they first induced 


star formation at the centres of their host galaxies about 8 Gyr ago. Thus, disc-like 
bulges with stars as old as 8 Gyr are perfectly possible. 

For the disc to be replenished with gas so that the bar can push this gas to the 
centre and produce a new central burst of star formation and rebuilding of the disc¬ 
like bulge, it has to fall into the disc plane from a direction not parallel to the galaxy 
disc, and inside the bar corotation radius. Otherwise, this gas will be pushed out¬ 


wards by the bar or accumulate at the corotation (see Bournaud and Combes 2002). 


Evidence for gas infall from directions not parallel to the galaxy disc has recently 
been presented by Bouche et al ( 2013} , but how often it occurs is still unknown, as 
is whether the gas reaches the disc within corotation. 

Thus, although Ellison et alj( [2011} find that, statistically, barred galaxies present 
ongoing, central star formation more often than unbarred galaxies, there is still a 
significant fraction of barred galaxies with star formation rates comparable to those 
in unbarred galaxies. In addition, Coelho and Gadotti (2011) find that the younger 
bulges found in barred galaxies have a mean stellar age of a few Gyr. This is in 
contrast to unbarred galaxies, which show on average older mean stellar ages (see 
Fig. |T0]>. This means that replenishing the disc with gas inside the corotation radius 
is a phenomenon that does not occur very often. Otherwise, very young stellar pop¬ 
ulations should be more conspicuous at the centres of barred galaxies. As discussed 
above, Bensby et al ( 2013b} ) find a stellar component in the Milky Way bulge with 
ages less than 5 Gyr. Although this means that the mean stellar age for the Milky 
Way bulge as a whole is above this value, this younger component has a mean stel¬ 
lar age thus in very good agreement with the mean stellar ages of the young bulges 
in other barred galaxies. The bottom line is that stars originating from gas infall to 
the centre through disc instabilities do not necessarily have to be extremely young 
now. A fraction of these younger stars could be elevated out of the disc plane and 
populate the box/peanut, but most of these stars are expected to be at or near the 
large scale disc plane. 


- Gas outside corotation receives angular momentum from the bar, and other factors govern the gas 
infall rate at these outer radii, such as dynamical effects induced by spiral arms and the dissipative 
nature of the ISM. At t hese d istances the infall of gas no longer occurs in a dynamical time scale. 

3 Some simulations iKraljic et al 2012) have reported the early formation of bars, at redshifts 
above 1. However, these bars are short-lived. In these simulations, bars formed at < 1 generally 
persist down to z = 0. 
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Fig. 10 Relative distributions of bulge mean stellar ages for barred and unbarred galaxies in bins 
of bulge stellar mass, as indicated. For massive bulges, the distribution is bimodal only for barred 
galaxies, consistent with the picture in which secular evolution processes build disc-like bulges. 
However, note that the mean stellar ages of such bulges can be as high as a few Gyr. [Adapted from 
\Coelho and Gadotti\2011\ \ 


Some studies using both simulations and observations have suggested that some 
star formation may occur along the bar as soon as the bar forms and induces shocks 


in the gas content in the interface between the bar and the disc (see Athanassoula 


1992 1 Phillips 1996 Sheth et al 2002!. However, they also indicate that, due to 


shearing in the gas clouds when they start falling towards the centre, this is inter¬ 
rupted shortly afterwards, and from then on, star formation is limited to components 
at the bar ends and the disc centre. Discs stars trapped by the bar mostly do not leave 
the bar, keeping their elongated orbits. As it evolves, however, the bar can capture 
stars from the disc formed more recently (see an example of such capture in the 
right panels of Fig. [8j, and thus a fraction of young stars can be present in the bar, as 
long as there is ongoing star formation in the disc within the radius at which the bar 
ends, and the bar can grow stronger and keep capturing disc stars. Stars in bars are 
thus predominantly old, and therefore box/peanuts, being just part of bars, should 
as well be populated mostly by old stars, with some younger component. 

For a galaxy as massive as our own, the bar is expected to form at redshifts 
close to 1 (see |Sheth et~al) |2012j i. The ages of the stellar populations seen at the 
Galactic bulge are thus consistent with the picture of it being built purely from the 
bar instability in the disc, i.e. the Galactic bulge can well be just a box/peanut plus 
a disc-like bulge, as far as the ages of its stellar populations is concerned. 

In trying to assess how the bulge of the Galaxy has formed, the chemical content 
of its stellar populations are better consider closely with their kinematical properties. 
Bars and their box/peanuts are expected to rotate cylindrically, i.e. the mean stellar 
velocity is independent of the height above the plane of the disc - as a rigid body. 
Variations in this pattern are usually attributed to the presence of other structural 
components with different kinematical properties (see e.g. Williams et al 2011 and 
references therein). As discussed above, one stellar population in the Bulge can be 
described as having high metallicity, low content of alpha-elements, and kinematics 
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consistent with the eccentric stellar orbits in bars in cylindrical rotation. The alpha- 
element content thus indicates that the box/peanut formed after the thick disc. In 
external galaxies, recent evidence suggests that thick discs form early, during the 
short initial formation stages of the galaxy, which qualitatively agrees with the pic¬ 
ture for the Milky Way (Comeron et al] 201l||2014| i. Nevertheless, although most 
of the thick disc stars thus form in situ, a significant fraction of stars in the thick 
disc may come from the accretion of satellites, and another fraction (likely smaller), 
may consist of heated up stars from the thin disc. Evidently, this complicates sub¬ 
stantially the interpretation of observations. 

Likewise, the observations from the ARGUS survey indicating variations in 
metallicity within the box/peanut (see Sect. 4.1), are suggestive of a complex stellar 
population content in the box/peanut itself. This can be a result of similar processes 
that complicate just as well the stellar population content of the thick disc, as de¬ 
scribed in the previous paragraph. As discussed in Sect. 4, different populations or 
population gradients in the disc from which the box/peanut forms can produce a 
similar result. However, it can also be the result of more than one buckling event 
forming the box/peanut. In the simulations of Martinez-Valpuesta et alj( ]2006) a bar 
goes through a first buckling event about 2 Gyr after the beginning of the simulation, 
and this event is fast (< 1 Gyr). However, a second, powerful buckling event occurs 
around 5 Gyr later and lasts for about 3 Gyr. How this would affect the chemical 
content of stars seen today in the box/peanut depends on how the stellar population 
content and kinematical properties of the bar and box/peanut vary during these peri¬ 
ods. But it is clear that if the bar of the Milky Way has gone through such recurrent 
buckling, the presence of multiple populations in the box/peanut is not surprising. 
Add this to the complex composition of the thick disc and one sees how complicated 
the stellar populations can be away from the disc plane. 


7.2.2 Bulges formed via violent processes 


The classical picture of bulges in disc galaxies is that of mini-ellipticals: massive, 
smooth and extended spheroids with dense centres, with old stellar populations 
showing alpha-enhanced chemical composition and relatively low rotational sup¬ 
port (as compared to discs). A natural formation scenario for these structural com¬ 
ponents would be that of a monolithic collapse (Eggen et al, 19621, in which a single 
gas cloud collapses in short time scales (< 10 8 yr), producing a violent burst of star 
formation that originates the stellar halo and the bulge. While this scenario might 
explain the formation of the first spheroids, it faces many difficulties. It does not re¬ 
produce for instance the heterogenous distributions of stellar ages and metallicities 
observed in bulges (e.g. McWilliam and Rich 1994[ Wyse et al 1997 a monolithic 
collapse implies a more homogeneous population) and regions of ongoing star for¬ 
mation (e.g. Carollo et al 1997| . 

It should be noted, however, that the monolithic collapse scenario was formulated 
within a perspective that does not include box/peanuts, and thus should not be com¬ 
pared against the properties of such observed structures. On the other hand, modern 
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dissipative collapse models similar to the pioneer model of |Eggen et al| ( jT%2[ ), that 
however include as well cosmological ingredients, generate bulges with more real¬ 
istic properties. Such key ingredients include in particular a long time scale history 
of accretion of dark matter haloes into the central halo, with the associated evolution 


of angular momentum and star formation episodes (see e.g. Samland and Gerhard 
[2003] IQbreja etaH|20l3> . 

A natural picture within ACDM cosmology is that of merger-built bulges. Brooks 
& Christensen (this volume) review this picture. Another scenario to explain the 
formation of classical bulges is the coalescence of clumps in discs at high redshifts. 
This is also reviewed in this volume by Bournaud. 

The evidence from observations of the Galactic bulge as reviewed above, how¬ 
ever, give little support to the presence of a massive classical bulge. The observed 
box/peanut and its cylindrical rotation cannot be originated in violent scenarios. 
On the other hand, the possibility of a small classical bulge embedded within the 
box/peanut is not yet ruled out. Such composite bulges are discussed in the next 
section. 

Nevertheless, we have seen above that there is a component in the Galactic bulge 
with low metallicity and an alpha-element content consistent with it being formed 
concomitantly with the thick disc, i.e. before the box/peanut. In addition, the mor¬ 
phological properties of this component seem to point out a more spherically dis¬ 
tributed structure (see Fig.jTj). This is revealed by RR Lyrae stars and are properties 
that are shared by classical bulges. This component has a spatial extent similar to 
that of the box/peanut, but it is not revealed by images such as those from COBE. 
However, a number of early-type disc galaxies with massive classical bulges shows 
bars, which probably went through a buckling process that originated a box/peanut. 
So it is perfectly plausible to have a classical bulge and a box/peanut coexisting in 
the same galaxy. A possible example is our own massive neighbor, M31. Athanas-) 
|soula and Beaton| ([2006 ) have shown that this galaxy has a bar and a box/peanut, 
and there is evidence that it also hosts a classical bulge (e.g. Courteau et al| 201 1[ 
and references therein). 2D decompositions in Gadotti & Erwin (in preparation) 
show that the classical bulge has a similar extent as the box/peanut. The morphol¬ 
ogy of the M31 bulge (see Fig. 2 in Athanassoula and Beaton f2006 1 is similar to 
that seen in the COBE/DIRBE image for the Milky Way, although the vertices of 
the box/peanut are more clearly recognized in the Galaxy (perhaps due to projec¬ 
tion effects in M31). Nevertheless the same 2D decompositions of Gadotti & Erwin 
reveal the X shape outstandingly (see Gadotti 2012 1 . Another critical issue is the 
understanding of how bright/massive is the component revealed by RR Lyrae stars 
in the Galaxy, and how does this compare to other classical bulges. It clearly cannot 
be large enough as to mask the vertices of the box/peanut revealed by COBE. 


7.2.3 Composite bulges 

In the previous subsections, we explored the possibility that the Milky Way has a 
disc-like bulge, apart from its box/peanut. We also remarked about the possibility 
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of a classical bulge. Here we will briefly summarise recent work on the presence 
of such composite bulges in external galaxies. It is not hard to contemplate the 
possibility of such composite systems. A pure disc galaxy can form at high redshifts, 
say z ~ 3, and acquire a classical bulge, be it through minor mergers, accretion 
events or the coalescence of clumps of stars and gas in the disc. At z ~ 1 the same 
disc - now hosting a classical bulge at its centre - might become unstable to the 
formation of a bar, and develop one. Quickly this bar pushes gas to the central 
regions of the disc, originating a disc-like bulge. Give it a couple of Gyr and the 
box/peanut is formed. We thus end up with a galaxy containing a classical bulge, a 
disc-like bulge, and a box/peanut. 

|Gadotti| ( |20091 > has shown evidence that 34 per cent of his disc galaxies host¬ 
ing classical bulges are galaxies possessing bulges with structural properties typical 
of classical bulges, but with an intensity of star formation activity characteristic of 
disc-like bulges. While one cannot rule out the possibility that some classical bulges 
may present ongoing star formation, it is also plausible that many of these galaxies 
actually host composite bulges. Those would be composed by an extended classical 
bulge with an embedded disc-like bulge. Because the classical bulge dominates the 
disc-like bulge in terms of mass, the composite bulge shows structural properties 
of classical bulges. However, using the D„(4000) spectral index, Gadotti (2009) was 
able to realize the intense star formation in the bulge region. Such star-forming activ¬ 
ity, according to this interpretation, occurs at the embedded disc-like bulge. Further 


evidence for this and other types of composite bulges is presented in Mendez-Abreu 
|eTaI|(|20T4j). 


Kormendy and Barentine (|20 10 1 report the existence of a small disc-like bulge 


embedded in the box/peanut of NGC 4565. Nowak et~al| ( 2010( ) argue that NGC 3368 
and NGC 3489 actually have an embedded classical bulge within component(s) built 
via disc instabilities. Finally, very recently, Erwin et ~aT| (2014 > have shown further 
evidence of such components that appear to be embedded classical bulges. 

It will thus be no surprise if the Galactic bulge is a composite bulge. 


8 Concluding remarks 

In this review we have presented a summary of current progress towards character¬ 
ising the properties of the Milky Way bulge. In recent years, the spectroscopic and 
photometric surveys of the Bulge have provided us with the necessary tools to build 
a bridge connecting the detailed stellar population properties with a global view 
of the Galactic bulge. As a consequence, it is now becoming possible to discuss 
the Bulge properties as seen from an extragalactic perspective. Such a comparison, 
powered by the increasing number of models to which observations can now be di¬ 
rectly compared, is the only way in which we can set the history of events that led 
to the properties of the Bulge we see today. 

The era of surveys looking towards the Galactic bulge was born not only from 
our intrinsic desire to explore, but also as a response to the increasing complexity in 
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Bulge properties revealed with previous individual observations. The need to further 
map its morphology, to better constrain the spatial variations of its stellar popula¬ 
tions and subsequently connect it with their kinematics required the larger spatial 
coverage offered by such large-scale observations. 

While the dominant B/P nature of the Milky Way bulge has now been well estab¬ 
lished, it remains to be understood if the observed stellar population properties relate 
solely to the same structure or if they each have a different origin. Another way to 
phrase this would be to ask the following question: do the metal-poor, a-element 
enhanced, old bulge stars belong to a different structure than the B/P, which was 
formed somehow independently to the buckling instability process of the bar (i.e. 
as a classical bulge)? Currently, while we have important evidence from the con¬ 
nections between kinematics and metallicities of Bulge stars as well as their spatial 
distribution, we can only suspect about the presence of different components. How¬ 
ever, creating the link between these components and the specific formation sce¬ 
narios should be done with extreme caution, as a number of processes could have 
played a specific role at different stages of the assembly of the Galaxy. For example, 
the coalescence of disc clumps and the accretion of gas could have formed a thin 
disc, a thick disc, and even a spheroid in the centre during the early stages of for¬ 
mation of the Milky Way. The merging history could have also contributed to this 
assembly, which will depend on gas content, mass ratio and orbital parameters of 
the mergers, until the formation of a bar and the onset of the buckling instability 
took place to shape the dominant central component we see today - the B/P Bulge. 
Only by carrying out an extensive comparison of all the observational properties of 
Bulge stars with models and external galaxies, can we constrain the importance of 
all these events during the formation history of the Bulge. 

Acknowledgements We are grateful to an anonymous referee for many useful comments. We 
warmly thank Istvan Dekany for kindly providing us the table of individual distances to the RR- 
lyrae that we used to produce Figure 1. 
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